Realistic two-nucleon potentials for the relativistic two-nucleon Schroedinger equation 
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The potentials V(v) in the nonrelativistic ( relativistic) nucleon-nucleon (NN) Schrodinger equa- 
tion are related by a quadratic equation. That equation is numerically solved, thus providing phase 
equivalent v - potentials related for instance to the high precision NN potentials, which are adjusted 
to NN phase shift and mixing parameters in a nonrelativistic Schrodinger equation. The relativistic 
NN potentials embedded in a three-nucleon (3N) system for total NN momenta different from zero 
are also constructed in a numerically precise manner. They enter into the relativistic interacting 3N 
mass operator, which is needed for relativistic 3N calculations for bound and scattering states. 

PACS numbers: 21.45+v 



I. INTRODUCTION 

Traditionally the (semi) phenomenological high pre- 
cision two-nucleon (NN)potentials AV18 [H, CD Bonn 
and Nijm 1,11 H go together with the nonrela- 

f 2 

tivistic operator for the kinetic energy — in the NN 
cm. system. Nevertheless, as is well known, this non- 
relativistic Schrodinger equation 



(_ + v)tf 



k 2 



(1) 



can be related to an underlying relativistic NN 
Schrodinger equation 



{2\J m 2 + k 2 + v)<S> = 2 



fc 2 $ 



(2) 



by a simple algebraic step[3],[H|. Applying (2v m 2 + k 2 + 
v) to © from the left one obtains 

(4(m 2 + P) + 2uo{k)v + 2vuj(k) + P)$ = 4(m 2 + k%)<S>(3) 
which can be identically rewritten into {]]) if one defines 



V 



4m 



{2u{k)v + 2vuj(k) + v 2 



(4) 



with Lo(k) = v w 2 + k 2 . ( We use k in order to dis- 
tinguish the momentum operator from the number k). 
Therefore adjusting V in ([I]) to the NN phase shift and 
mixing parameters from a phase shift analysis and relat- 
ing the cm. momentum fco to the Lorentz invariant lab 
energy T iab via 



kn = 



mTiab 



(5) 
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( a relation identically valid for relativistic and nonrel- 
ativistic kinematics) one has in fact solved a relativistic 
equation. We also see that \& equals <E> .The question re- 
mains, what is v given V? The formal solution of that 
quadratic equation ((4]) is 



v = \jAmV + Auj(k) 2 - 2u>(k) 



(6) 



Why is v of interest? If one turns to the 3N system 
and would like to investigate relativistic effects @|,[3] the 
knowledge of v is very useful. It defines together with the 
relativistic kinetic energy the interacting NN mass opera- 
tor, which is a key ingredient for building the interactive 
3N mass operator [8(. Therefore we focus in this paper 
on the determination of v related to the high precision 
NN potentials via (gj) or @. 

In a potential v has been determined directly 

fitting to NN phase shifts. Thereby the Urbana V14 
potential has been readjusted achieving a fair fit ( though 
not of the quality of the high precision potentials). In 
lllll a momentum scale transformation 
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has been introduced which provides an analytical rela- 
tion between V and v and guarantees that the S-matrix 
related to |T]) at cm. momentum k equals the S-matrix re- 
lated to @ at cm. momentum q. In other words the rel- 
ativistic and nonrelativistic S-matrices agree at the same 
energy. This, however, is misleading since the equality of 
the two S-matrices should hold at the same cm. momenta 
12]. A better, though still approximate approach to re- 
late V and v has been given in Q. 

On the other hand there is the possibility to add an 
interaction to the square of the free NN mass operator 
h 2 = (2bj{k)) 2 +AmV ' . TheniT = |^-m has exactly the 
same form as ((T]) with V ' = V provided we identify h 2 
with the mass operator of the interacting NN system [i| , 
. This is of course also obvious from the relations ([T]) - 
( 2|) . The construction of the relativistic 3N Hamiltonian 



2 



requires, however, the 3N mass operator h rather than 
its square. Therefore our aim here is to solve ([U and © 
exactly for v. This is outlaid in section II. The validity 
of the resulting v is verified by demonstrating that it 
provides exactly the same phase shift parameters using 
([2|) as the underlying V using (TTJ). 

Next we regard the 3N mass operator where NN cm. 
forces enter in the form [8| 



)2 +p 2 



(2w(fc)) 2 + P 2 



(8) 



with p the total NN momentum. The p-dependence arises 
since in a 3-body system the NN subsystems are not at 
rest. In section III we propose a simple manner to deter- 
mine v p in a numerical precise way. This opens now the 
door to use v's which are equivalent to the underlying 
high precision potentials in a relativistic context in 3N 
bound and scattering problems. We summarize in sec- 
tion IV. A technical derivation is given in the Appendix. 



II. THE POTENTIAL V 

The determination of v using Eq. © can be achieved 
by a spectral decomposition. One can proceed in close 
analogy to the representation derived in [l3j for v p given 
in Eq. ([5]). We regard a specific partial wave state ( or 
coupled ones) with given orbital angular momentum(a) 
, total spin and total angular momentum. For the sake 
of a simpler notation we will not show these quantum 
numbers explicitely. Using the completeness relation of 
bound and scattering states for the potential V one ob- 
tains 



< k\v\k ' >=< k\^ b > M b < * b |fc ' > 
„ „J(k-k') 



+ / dk'\k "f < k\V'l > 2^k" 2 + m 2 < &{.\k ' > 
Jo 

(9) 

where ^b(k) is the nonrelativistic deuteron wave function 
of fTJ) and M}, the mass of the deuteron. Here we intro- 
duced a different definition of the binding energy, namely 
an implicit one: 



M b = \/An\ 2 + 4e b 'm = 2 



m + eb - 



(b 



4m 



(10) 



In lowest order it agrees with the usual one Mb — 2m +e/,. 
This new definition of the binding energy has in addition 
the feature that it can naturally be written as eb ' = — ^7 
in agreement with the form of the energy eigenvalue of 
([1]) at the bound state pole k = in of the S-matrix. 

The expression (|9]) can be identically rewritten into the 
form 

(k\v\k') = ^b{k)M b ^ b {k') 



-5{2uj(k)$l[T(k',k;—)} 
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T (k,k";—)T*(k',k";—)}. (11) 
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where T is the standard NN T-matrix related to V via 
the nonrelativistic Lippmann Schwinger equation. The 
derivation of that form is defered to the Appendix. A 
numerical implementation has not yet been performed, 
but we expect no problem. 

A second more simple way is to directly solve the 
quadratic operator equation (|4|). In momentum space 
it reads 

2m < k\V\k ' >= {uj{k)+uo{k ')) < k\v\k ' > 

+ - / dk "k " 2 < k\v\k " >< k "\v\k ' > (12) 
2 Jo 



< k\v\k ' > 

' ' dk "k " 2 < k\v\k " >< k "\v\k ' > 



2(Lu(k)+uj(k ')) Jo 
2m < k\V\k ' > 

Uj(k)+L0(k ') 



(13) 



We verified numerically that for all the realistic high 
precision potentials AV18, CD Bonn, and Nijm I, II the 
following very simple iterative scheme works 



< fc|v|fc ' > (0) = 



2m < k\V\k ' > 

Uj{k)+Uj(k ') ' 



(14) 



< fc|v|A; ' >(" +1 ) = 



2(w(fc) +cu(k ')) 



— x {4m < k\V\k ' > 



dk "k " 2 < k\v\k " >W< k "\v\k ' >W} (15) 



with n = 0, 1, ■ • ■ . For certain partial waves the iteration 
to converge requires an additional step, namely 

v(n+1) ^defme {av{n+1} + bv (n) )/{a + fc) ^ 

where the constants a and b are typically 1. 

We display in Table U an example documenting the 
convergence. In Figs [T][3] we show for an example the 
original nonrelativistic potential V(k, k '), the result- 
ing relativistic potential v(k, k ') and the difference 
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FIG. 1: The nonrelativistic potential V(k,k') AV18 in the 
state So- 




FIG. 2: The relativistic potential v(k, k') related to AV18 in 
the state 1 So- 



V(k,k ') — v(k,k '). We see in that example that 
V(k,k') > v(k,k'). 

To further characterize the difference between V(k, k ') 
and v(k,k ') one can regard the asymptotic behavior 

, /x const . , 

lim v {k,k') = — — (17) 

fe^oo,fc 'fixed k n 

against 

lim V(k,k') = C -^ (18) 

k— >oo,/c 'fixed K 

As examples , for AV18 we find n = 3, N = 2, whereas 
CD Bonn delivers n = 4.5, N = 3.5. The value of N for 
CD Bonn can easily be understood [lj] : N = 2+2+l/2- 
1 with the two 2's resulting from the meson propagator 
and the choice of the strong form factor, the 1/2 arising 
from transforming the Blankenbeclar-Sugar equation to 
the nonrelativistic Lippmann Schwinger equation and the 
(— 1) from the two Dirac spinors. In both cases, AV18 
and CD Bonn, n is larger than N by one unit, which is 
suggested by Eq. (fl~2|) . 

Having v at our disposal one can solve the relativistic 
Lippmann Schwinger equation 




4 \_ "W- _/>^ 
k[fm'] 6 8<^^2 4 k'tfin 1 ] 



FIG. 3: The difference between the relativistic and nonrela- 
tivistic potentials: V(k,k') — v(k,k') in the state So- 



TABLE I: Convergence of f/ n) to the iteration in Eq. (JT5J. 
We choose the coupled partial waves (^i^Di) of the Ar gonne 
V18 potential [l]]. The momenta k and k' are 1.0 fm _1 and 
the potential unit is [fm 2 ]. 



n 


» w ( 3 Si- s Si) 


w (n) ( 3 5i- 3 Di) 


v (n) ( 3 _Di- 3 Di) 





0.084232 


0.044709 


0.016853 


1 


0.067716 


0.044628 


0.016785 


2 


0.059933 


0.044597 


0.016744 


3 


0.056135 


0.044587 


0.016719 


4 


0.054234 


0.044585 


0.016705 


5 


0.053259 


0.044587 


0.016696 


6 


0.052749 


0.044589 


0.016691 


10 


0.052194 


0.044595 


0.016684 


20 


0.052126 


0.044597 


0.016684 


30 


0.052126 


0.044597 


0.016684 



t(k,k ') = v(k,k') 

+ / dk"k" 2 v(k,k") 
Jo 

MfcO-M*")*^ *'' (19) 

for the half shell t-matrix. It is related to the S-matrix 
via 

s(k) = e 2iSAk) = 1 - inkuj(k)t(k, k). (20) 

The corresponding nonrelativistic relation is 

S(k) = e 2iSnr W = 1 - iirkmT{k, k). (21) 

where T(fc, k ') obeys the standard nonrelativistic Lipp- 
mann Schwinger equation. We illustrate in Table [II] for 
the ( 3 Si - 3 Di ) partial wave states the perfect agreement 
of the corresponding phase shift and mixing parameters. 
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Comparison of the phase shift and the mixing 



TABLE II: 
parameters 

for the coupled partial waves ( 3 Si- 3 .Di) of the Ar gonne V18 
potential^. The second column points to the relativistic 
(Rel.) or nonrelativistic (Nonrel.) calculations. The unit 
of the phases are in degrees. 



b . [MeV] 


Nonrel./Rel. 


S( 3 S 1 ) 


H 3 Di) 


£ 


1.0 


Nonrel. 


147.62 


-0.0050743 


0.10303 


1.0 


Rel. 


147.62 


-0.0050744 


0.10304 


10.0 


Nonrel. 


102.71 


-0.66593 


1.1267 


10.0 


Rel. 


102.71 


-0.66593 


1.1267 


50.0 


Nonrel. 


62.929 


-6.3189 


2.0853 


50.0 


Rel. 


62.929 


-6.3189 


2.0854 


100.0 


Nonrel. 


43.531 


-12.093 


2.4899 


100.0 


Rel. 


43.531 


-12.093 


2.4899 


350.0 


Nonrel. 


2.6451 


-26.719 


4.9223 


350.0 


Rel. 


2.6452 


-26.719 


4.9222 



III. THE CONSTRUCTION OF THE NN FORCE 
EMBEDDED IN THE 3N SYSTEM 



If one defines 2ui p (k) = y A(uo(k)) 2 + p 2 then using (j4|) 
Eq. ((HI) can be written 



(v p + 2uj p {k)) 2 = (2a; (A;) + v) 2 + p 2 
4(w p (fc)) 2 + AmV 



or 



(v p ) 2 + 2v p ujp(k) + 2id p (k)vp = AmV 

Between free states it yields 

v p (k,k '){ojp{k) +ui p (k ')) 
= 2mV(k,k ') 

1 p OO 

' dk"(k") 2 V p (k,k")Vp(k" 1 k') 



(22) 



(23) 



(24) 



This has the same structure as (fl"4"]) replacing oj(k) by 
u)p(k). The iteration procedure described in section II 
works equally well. We compare in Fig H «(£;,£; ') and 
v p (k,k ')■ We see a weakening of v p against v. This is a 
fact known from previous calculations [l~3T ] and from the 
approximate ( but very useful) expression 



v p (k,k ') =v{k,k ')(! 



8u(k)(j(k ') 



77) 



(25) 



This is demonstrated in Fig[4l where also an even simpler 
approximation 



v p (k, k ') — v{k, k )(1 



P 

8m' 







(26) 



is shown. That latter approximate version, however, is 
somewhat worse and is not recommended. 
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FIG. 4: Comparison v(k,k ') against v p (k,k ') for AV18 in 
the state 1 S'o for fixed k' — 1.0 fm -1 . The solid curve shows 
the relativistic potential v (p — 0). The other curves show v v 
for p = 4.0fm _1 . The long dashed curve is the exact one, the 
short dashed and dotted curves show the approximations for 
v p given in Eqs. (|25[) and (|26[1 . respectively. 



IV. SUMMARY 

Relativistic calculations in the instant form of dynam- 
ics proposed in [l5[ requires an interacting 3N mass op- 
erator. As has been shown in 0] and used in [fl] . fj\ , [l3] . 
[IBli [I3|' an( i [9 the NN potential in a moving frame 
Vp enters in the form given in Eq. (jHJl , where the NN force 
in the NN cm. system v enters into the relativistic NN 
Schrodinger equation . We showed that the quadratic 
operator relation ^ for v can be solved directly in an 
iterative manner and this very precisely. This has been 
documented by evaluating NN phase shift and mixing pa- 
rameters using the standard nonrelatvistic Schrodinger 
equation ([TJ) and the relativistic NN Schrodinger equa- 
tion @. This opens the way to get w's related to any 
NN potential V adjusted in a nonrelativistic frame work 
like the high precision NN potentials. By the same itera- 
tive procedure also v p can be gained. It turned out that 
v is smaller in magnitude than V and v p is smaller than 
v. Applications to the 3N bound and scattering states 
are planned. 

The numerical calculations has been performed on the 
IBM Regatta p690+ of the NIC in Julich, Germany. 



V. APPENDIX 

Derivation of Eq. (fTTj) 
We start from ^ and use the well known decomposition 

, , , S(k-k') T(k,k';^-) , s 



L.2 

It - — 



to arrive at 



< k\v\k ' >= V b (k)M b V b (k ') 
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13 

- ie 



2co(k') 



T(h k"- 

dk'\k "f ir : m J 



T*(k',k";!£) 

k" 2 „•, k' 2 



2to{k") 



(28) 



The integral requires some care and we keep the limit- 
ing processes for the two scattering states separately by 
putting 

1 1 



k" 2 -k 2 +ie k" 2 - k' 2 - ie 
1 1 

k" 2 -k 2 + iei k" 2 - k' 2 - ie 2 
1 1 



k" 2 - k 2 + iei 
1 



k" 2 - k' z - ie 2 



(29) 



k 2 - k' 2 - t(e 1 + e 2 )' 
This allows us to perform one limit firstly with the result 
f 1 



k" 2 -k 2 +ie k" 2 - k' 2 - ie 

^ ■ r,,„2 ,2 



k" 2 - k 2 



nrS(k" z - k 2 ) 



1 



- ( ? , + ™S(k" 2 - k' 2 )) \ 1 

\k" 2 -k' 2 X 7 k 2 -k' 2 -ie 1 

Thus we get for some well behaved function f(k") 

f(k") 



(30) 



dk "(k "Y 



i 



(k" 2 - k 2 + iei){k" 2 - k'- ie 2 ) 



lim - 

e ^ +0 f.2 - f~' - it 



and therefore 



dk"(k ") z f{k")8{k" z - k 2 ) 



k 



112 



(33) 



= 2m 2 dk"(k") 2 cj(k")T{k,k"; ) 

Jo yn m 

T *(k' 1 k"-—)8{k" 2 -k 2 ) 
m 

k 2 k 2 
= m 2 kuj(k)T(k,k;—)T*(k',k;—) 
m m 



This is part of the unitary relation 



T(k,k'-X)-T*(k',k-X) 
m m 
a 2 

= 2i%T(k,k'; — ) 
m 

Q 2 

= 2i%T(k',k;—) 
m 

r°° a 2 a 2 

- 2iirm / dk"(k ") 2 T(k, k"; ?-)T*(k', k"; ?-)S(q 2 - k" 2 ) 
Jo mm 

2 2 

= —inmqT(k 1 q; — )T*(k',q; — ) 
m m 



(34) 



( Note we used the symmetry of the off- the- energy- shell 
T- matrix). Consequently 



dk"k " 2 f(k ")8{k" 2 - k ) 



2 -^Z[T(k,k'-X) 
7r m 



(35) 



and 



V 



dk"(k 



fell2 _ £,2 



-V 



f(k") 



dk"{k "Y 



lim 



1 



—in lim , x 

e ^+o fc2 _ k a _ ie 

POO 

( / dk"(k ") 2 f(k")s(k" 2 - k 2 ) 

Jo 

+ / dk"(k" ff{k")5{k" 2 -k' 2 )) 



;+o fc2 _ _ ie 
( / dk"(k "ff{k")5{k" 2 - k 2 ) 

dk"(k ") 2 f{k")5{k" 2 -k' 2 )) 
2mi 

hm 



(31) 



e— +o k 2 - k' 2 - ie 

f k 2 k' 2 \ 

w(k)$S[T(k, k'; —)] + <j{k')S[T(k k; —)] 36) 
m ml 



The principal value prescription is denoted as "V J" . In 
our case 

i m u/2 Combined with Eq. (|28[) certain terms cancel and one 

f(k") = 2m 2 oj(k")T(k 1 k"; )T*(k',k"; ) (32) arrives at Eq.jm 

m m 
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